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Abstract 

We derive Feynman rules for the interactions of a single gravitino with (s) quarks and glu- 
ons/gluinos from an effective supergravity Lagrangian in non-derivative form and use them to 
calculate the hadroproduction cross sections and decay widths of single gravitinos. We confirm 
the results obtained previously with a derivative Lagrangian as well as those obtained with the 
non-derivative Lagrangian in the high-energy limit and elaborate on the connection between gauge 
independence and the presence of quartic vertices. We perform extensive numerical studies of 
branching ratios, total cross sections, and transverse-momentum spectra at the Tevatron and the 
LHC. From the latest CDF monojet cross section limit, we derive a new and robust exclusion con- 
tour in the gravitino-squark/gluino mass plane, implying that gravitinos with masses below 2- 10 -5 
to 1 • 10 -5 eV are excluded for squark/gluino- masses below 200 and 500 GeV, respectively. These 
limits are complementary to the one obtained by the CDF collaboration, 1.1 • 10~ 5 eV, under the 
assumption of infinitely heavy squarks and gluinos. For the LHC, we conclude that SUSY scenarios 
with light gravitinos will lead to a striking monojet signal very quickly after its startup. 
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I. INTRODUCTION 



Together with the possible existence of extra spatial dimensions, supersymmetry (SUSY) 
remains the prime candidate for physics beyond the Standard Model (SM). Among the 
many undisputed theoretical advantages of the Minimal Supersymmetric SM (MSSM), the 
intimate connection of this new space-time symmetry with electroweak symmetry breaking 
is of particular importance. The search for SM or MSSM Higgs bosons as well as for spin-0 
and spin- 1/2 partners of the SM fermions and gauge bosons are therefore often considered 
to be the most important tasks for present and future collider experiments. 

For many years, the focus has been on minimal supergravity (mSUGRA) models, in which 
SUSY is broken by gravitational interactions and the lightest SUSY particle (LSP) is the 
photino or, more generally, the lightest of four neutralinos, Xv Only around 1980 it was 
discovered that the SUSY partner of the spin-2 graviton, the massless spin-3/2 gravitino, 
does not necessarily couple to matter with gravitational strength only, but that its coupling 
can be enhanced to electroweak strength once SUSY is broken through the super-Higgs 
mechanism and the associated Goldstone fermion, the spin-1/2 goldstino, is absorbed to 
give the gravitino its mass and its longitudinal degrees of freedom, making it the LSP 

HQ. 

The electroweak strength of goldstino interactions with massless photons and photinos 
was then used to impose limits on the gravitino mass by comparing total theoretical cross 
sections for electron-positron colliders to experimental single-photon searches at PEP and 
PETRA, resulting in a first mass limit of mg > 2.3 • 10~ 6 eV [3[. Subsequently, the single- 
photon searches at LEP 1 and LEP 161 with cross section limits of 0.1 and 1 pb implied 
gravitino masses above 10~ 3 and 10 -5 eV for light neutralinos of mass below 50 and 100 
GeV, respectively 4j. These limits where, however, obtained without imposing missing or 
observed photon energy cuts on the theoretical cross section. 

In 1988, the CDF collaboration published a cross section limit of 100 pb for their monojet 
search at the Fermilab pp collider Tevatron [^j], which they used to impose bounds on the 
squark-gluino mass plane, but which could also be interpreted as the absence of a light 
(s) goldstino signal, yielding > 2.2 ■ 10~ 5 eV and m- g > 100 GeV [6J. This first hadron- 
collider analysis assumed, however, very heavy squarks of rriq > 500 GeV and was based 
on partonic subprocesses involving only gluons and gluinos, but no (s) quarks. The analysis 
was later re-applied to the 1996 CDF multijet cross section limit of 1.4 pb 7], yielding 
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m G — 3 ' 10 -4 eV and m- g > 200 GeV [8]. Predictions were also made for the LHC, albeit 
for an assumed center-of-mass energy of 16 TeV 9|. 

While the gluon-gluon initial state dominates indeed for the production of light final states 
at the LHC, it is well known that it is the quark- ant iquark luminosity that dominates at the 
Tevatron and that quark-gluon initiated QCD Compton processes contribute significantly 
for heavier final states at the LHC. A complete and robust study must therefore take into 
account 1) all partonic subprocesses leading to the production of single gravitinos, i.e. qq — > 
Gg, gg — » Gg, and qg — > Gq, 2) the subsequent decay of the squark/gluino into an observed 
jet and a second gravitino, i.e. g — > Gg and q — > Gq, 3) up-to-date collider energies, parton 
density functions (PDFs), values of Aqcd, and SUSY-breaking scenarios, 4) the experimental 
cuts on the jet and missing transverse energies, and 5) the most recent experimental cross 
section limits. 

In Sec. Ull we calculate analytical gravitino production cross sections and decay widths 
using an effective supergravity Lagrangian in four-component notation (see App. |A} and the 
Feynman rules for single gravitinos derived from it (see App. [B]). We discuss in some detail 
the gauge-independence of our results and its relation to the sign of interferences and the 
presence of quartic vertices. In Sec. IHIt we first present a concise review of gauge-mediated 
SUSY-breaking (GMSB) models, where gravitinos are naturally the lightest SUSY particles, 
and discuss their implementation in different benchmark slopes. We then establish the 
regions in which gluino/squark decays into gravitinos and jets dominate. Next, we present 
the various subprocess contributions to the total cross sections at the Tevatron and LHC 
and compute the jet and missing transverse-momentum spectra, taking into account the 
gluino/squark decays. Finally, we deduce a new limit on the gravitino mass from the latest 
CDF monojet search and discuss the signal size and missing-i^r trigger thresholds at the 
LHC. Our conclusions are presented in Sec. HVl The discussion of cosmological constraints 
on the gravitino mass is beyond the scope of this paper. For a recent analysis of Lyman-a 
forest and WMAP data, assuming a light gravitino as a warm dark matter candidate in 
GMSB models and yielding itlq < 16 eV, we refer the reader to [10'] and the references 
therein. 
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II. ANALYTICAL RESULTS 



In this Section, we present our analytical results for the hadroproduction cross sections 
of single gravitinos with gluinos and squarks (Sec. Ill Al) and the two-body decay widths of 
squarks and gluinos into gravitinos with quarks and gluons (Sec. Ill B[) . They have been 
obtained by using an effective supergravity Lagrangian in non-derivative form (see App. |A"1). 
from which the corresponding Feynman rules (see App. [Bj have been derived. 



A. Production 



In .R-parity conserving supersymmetry, single gravitinos can be produced in strong inter- 
actions in association with either gluinos or squarks. In addition, the associated production 
of gravitinos and gluinos proceeds through two competing initial states, i.e. quark- ant iquark 
or gluon-gluon scattering, while gravitinos and squarks can only be produced in quark-gluon 
scattering due to fermion number conservation. The differential cross sections 

- = -—W 2 (i) 

dt 2s 8ns 1 1 w 
depend in general on the SUSY particle masses, the usual Mandelstam variables s, t, and u, 
and their mass-subtracted counterparts, t^ g = t — m?^ and Uq s g = u — rn^ g . The gravitino 
mass ttLq will be neglected everywhere except in the coupling constants, so that t is integrated 
over the interval [— s + m~g ; 0]. 

We first consider the process initiated by quarks and anti-quarks, 

qq^G~g, (2) 

whose contributing Feynman diagrams are shown in Fig. [TJ The corresponding squared 
transition matrix element, averaged (summed) over initial (final) state spins and colors and 
summed over left- and right-handed squark exchanges, 

T 2 9 2 s C F 



-.2 /~i l~™2 ™4 ™4 



|M| 



qq->G~g 3 M 2 m 2 



— £ (2tu-m 2 Jt + u)) +-f^j + 4««s 
9 —^ (-2tu + m\(t + u) 



(3) 



is symmetric under the exchange of the t and u Mandelstam variables. The s-channel 
contribution is individually gauge- independent, and the t- and w-channel contributions are 




FIG. 1: Leading-order Feynman diagrams for the production of a gravitino in association with a 
gluino in quark-antiquark collisions. 



manifestly gauge-independent. For gluino pair production, an interference term between the 



t- and u-channel diagrams proportional to the squared gluino mass exists [ll|] , but this term 
vanishes for gravitino-gluino associated production linearly with the gravitino mass. 

The contributions from individual diagrams that we obtain differ, of course, from those 
presented in Eq. (4) of [l^ j. since our effective Feynman rules are proportional to the SUSY 
particle masses, but our total results agree. A related cross section has been computed with 
effective Feynman rules in Eq. (28) of [4] for the associated production of gravitinos and 
neutralinos at lepton colliders. After adjustment of masses, couplings, and color factors, it 
agrees with our result. In the limit of negligible squark- and gluino masses, where the t- 
and w-channel contributions vanish both due to the higher mass- dimension of the squark 

n 

coupling, our result agrees also with Tab. 1 in 13] when summed over left- and right-handed 



quarks. This limit is, however, only applicable in the high-energy context of cosmology 
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FIG. 2: Leading-order Feynman diagrams for the production of a gravitino in association with a 
gluino in gluon-gluon collisions. 

and not at current hadron colliders. 

Next, we compute the competing gluon-initiated process 



99 -> Gg, 



(4) 



whose contributing Feynman diagrams are shown in Fig. [2j The gauge-independent total 
squared matrix element 



|M| 



9s m § 



stu 



99->Gg QC F M 2 m 2 d S 2 t 2 gU 2 g 



tu (t 2 + u 2 ) - ml (t 3 + 6t 2 u + 6tu 2 + u 3 ) 



+2m 4 , (2t 2 + 7tu + 2m 2 ) - 5m? (t + u) 



(5) 



averaged (summed) over initial (final) state spins and colors, is again symmetric under inter- 
change of the final gravitino and gluino and consequently also of the Mandelstam variables 
t and u. Our result agrees with Eq. (6) in 12] and also with Tab. 1 in [13] in the limit of 
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FIG. 3: Leading-order Feynman diagrams for the production of a gravitino in association with a 
squark in quark-gluon collisions. 



small gluino mass. It also agrees with Eq. (5) in |8j, if its variables t and u are understood 
to be their mass-subtracted counterparts and if its integration variable z, whose definition 
is unfortunately missing, is assumed to be given by z = 1 + 2t/(s — m~) G [—1; 1]. 

Finally, we analyze the associated production of gravitinos and squarks, which is initiated 
by quark-gluon scattering, 

qg -> Gq h (6) 

and proceeds through the Feynman diagrams in Fig. [31 In this case, the squared matrix 
element for a squark of a given chirality i is 



M 



G L 



st% 



111 



2 4 



i(2u) 



(7) 



where the first term in the squared brackets denotes the gauge-independent sum of s- and 
t-channel contributions including their interference, while the u-channel contribution in the 
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FIG. 4: Leading-order Feynman diagrams for the decay of a squark into a gravitino and a quark 
(left) and a gluino into a gravitino and a gluon (right). 

second term is individually gauge- independent. The third term corresponds to the gauge- 
independent sum of s- and t-channel interferences with the w-channel. The same result is 
obtained for the charge-conjugated process, 



Q9 G Q* 



(8) 



This leads to a factor of two for pp colliders with a neutral initial state such as the Teva- 
tron, but not for pp colliders such as the LHC, where the parton densities are not charge- 
symmetric. Our result agrees with the one in 12j, which has been obtained with derivative 
couplings including a quark-gluon-gravitino-squark vertex contribution. Note that this quar- 
tic vertex is absent in the effective theory 14] . since it would spoil the gauge-independence. 



At high energies, only the w-channel contribution survives, so that our result agrees with 
the one in Tab. 1 of 131 ] . 



B. Decay 

Heavy squarks and gluinos may decay either directly or through cascades into the lightest 
SUSY particle, which we assume to be the gravitino. Direct decays, which dominate for light 
gravitinos [6|, proceed through the Feynman diagrams shown in Fig.Hl and the corresponding 
partial widths 

dr 1 1 



dt ^ m q,g 87rm| ^ 



\M\ (9) 



S 



are obtained from the squared transition matrix elements after integration of the Mandelstam 
variable t over the interval [— m|g + m 2 q g ; 0]. For a squark of a given chirality i, the squared 
matrix element 



7'2 



(m| - m|) 3 
4 



leads then to the partial width 
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Since the gluon mass is, of course, zero (m g = 0), the squared gluino decay matrix element, 
averaged (summed) over initial (final) spins is 

72 m% 



Or 



leading to the partial width 



Or 



15]. They agree, in particular, with the general result in Eq. 



These results are well-known 

(6.24) of [3], valid for the decay of any heavier SUSY particle into its Standard Model 
partner and a lighter gravitino. 

III. NUMERICAL RESULTS 

A. Gauge Mediated Super symmetry Breaking 

In Gauge Mediated Supersymmetry Breaking (GMSB) models, SUSY breaking occurs 
in a secluded sector at the scale (F), related to the gravitino mass by ttlq = (F)/(v / 3M), 
and is transmitted to the observable sector by a chiral superfield S and n q quark-like and ni 
lepton-like messenger fields HQ- The superfield S is a gauge-singlet, but its scalar and 
auxiliary components overlap with the gravitino and acquire vacuum expectation values 
(S) and (Fg). The messenger fields then acquire a mass M mcss ~ (S) through Yukawa 
couplings to the superfield S. They are given the same Standard Model gauge couplings to 
the observable fields as ordinary quarks and leptons, so that they can induce gaugino and 
sfermion masses through one- and two-loop self-energy diagrams, respectively. The lightest 



SUSY particle is always the gravitino, and it is for this reason, that we concentrate our 
numerical study of gravitino production at hadron colliders on GMSB scenarios. 

Besides M mess , n q , and ni, GMSB scenarios are determined by the ratio of Higgs vacuum 
expectation values, tan/3, the sign of the Higgs mass-parameter /i, and by the auxiliary 
vacuum expectation value, (Fg), which is related to the mass splitting of the messenger 
fields and realistically considerably smaller than both the squared mass scale of the messenger 
fields, (S) 2 , and the fundamental SUSY-breaking scale, (F). It is usually re-expressed in 
terms of an effective SUSY-breaking scale, A = (F S )/(S). 

A number of SUSY benchmark scenarios have been proposed in [3] in order to facilitate 
detailed comparisons between SUSY searches at different colliders and with different sig- 
nals/backgrounds. In particular, we show in Fig. [5] the average squark and physical gluino 
masses for the two GMSB scenarios proposed in 
effective SUSY-breaking scale A, with M mess /A = 2, tan/3 = 15, and /i > fixed. For 
SPS 7, the benchmark point (indicated by a vertical dashed line) is at A = 40 TeV and 
n q — ni — 3, leading to a stau (fi) next-to-lightest SUSY particle (NLSP), while for SPS 
8, the benchmark point is at A = 100 TeV and n q = m = 1, leading to a neutralino (x?) 
NLSP. The regions that have already been excluded by LEP2 and Tevatron searches for 
light neutralinos and charginos in GMSB scenarios lie to the left of the vertical dotted line 
[l9( . The physical masses in Fig. [5] have been obtained by imposing boundary conditions at 
the Grand Unification Theory (GUT) scale and evolving them to the electroweak symmetry 
breaking scale via renormalization group equations using the computer program SUSPECT 
20| | . Note that the mass hierarchy of gluinos and squarks at SPS 7, m g > Trig, is reversed 
at SPS 8, where m g < Trig. 



18], SPS 7 and 8, as a function of the 



B. Branching Ratios 

We are now in a position to determine the regions in SUSY parameter space, where the 
squarks and gluinos, that are produced in association with the gravitino at hadron colliders, 
decay dominantly into a two-body final state with a second gravitino and a quark or gluon, 
leading to an experimentally identifiable monojet signal with large missing transverse energy. 

To this end, we evaluate the decay widths calculated in Sec. Ill 131 in the GMSB 
scenarios discussed in Sec. IIII Al and compare them to the competing total decay width 
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FIG. 5: Average squark and physical gluino masses for the GMSB benchmark slopes SPS 7 (top) 
and SPS 8 (bottom) as a function of the effective SUSY-breaking scale A, with M mess /A = 2, 
tan/3 = 15, and \i > fixed. 
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Tmssm of gluinos and squarks into MSSM two-body final states up to one-loop level and 
those into three- and four-body final states at tree-level as implemented in the computer 
program SDECAY [2]]]. The resulting branching ratios 

BR = \ r (14) 

1 MSSM + 1 G 

of gluinos and squarks are shown in Figs. [6] and [7] as a function of the gravitino mass and 
of the effective SUSY-breaking scale A, defining the GMSB benchmark slopes SPS 7 (top) 
and SPS 8 (bottom). 

For both benchmark points, we observe that left- and right-handed, up- and down-type 
squarks decay dominantly with BR > 0.9 into gravitinos, if the gravitino mass does not 
exceed m G < 10 -4 eV. While the SUSY-breaking scale A must not lie significantly below 
the benchmark points of A = 40 and 100 TeV, respectively, these regions are already largely 



excluded 



jy LEP2 and Tevatron searches for light neutralinos and charginos in GMSB 



scenarios .19]. 



For gluinos, the conclusions are quite similar for SPS 7, but more optimistic for SPS 8 
with the decay into gravitinos dominating up to < 10~ 3 eV for all physical values of A. 
Above these limits, the decay chains depend essentially on the mass hierarchy of the SUSY 
spectrum with g qq and q — > X^q^ at SPS 7, whereas q — > gq and g — > Xi' QQ^ a ^ 
SPS 8, leading in general to more complicated multijet signals. Note that all of these decays 
are instantaneous with decay lengths around or below 1 fm, so that they occur close to the 
primary vertex and well inside any collider detector. 

While our numerical results for squark decays are new, gluino decays have been studied 

Siite some time ago in a simple SUSY scenario with a massless photino NLSP. In Fig. 1 of 
], the two-body decay g — > Gg has been compared with the tree- level three-body decay 
9 IQQi neglecting all other decay modes and for m g = 100 GeV and rriq = 500 GeV. The 
conclusion there was that BR > 0.9 up to m G < 10 -4 eV, which compares quite favorably 
with our result at SPS 8 and A = 20 TeV (see the lower part of Fig. [7]), where the masses 
trig = 200 GeV, = 250 GeV, and m^o = 14 GeV are of similar magnitude. For rriq = 1000 
GeV, the gravitino decay mode was found to dominate up to m G = 5 ■ 10~ 3 eV (see Fig. 1 
of js|). This compares again favorably with our result at the benchmark point SPS 8 (see 
lower part of Fig. EJ), where Trig = 1100 GeV. Related results for a massless photino NLSP 
and gluino masses between 200 and 750 GeV and squark masses between 500 and 2000 GeV 
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FIG. 6: Branching ratios of gluinos/squarks into gravitinos and gluons/quarks as a function of the 
gravitino mass with the other SUSY masses fixed to their SPS 7 (top) and SPS 8 (bottom) GMSB 
benchmark values. 
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can furthermore be found in Fig. 1 of 9]. 



C. Tevatron 

The total hadronic cross section for gravitino-gluino or gravitino-squark associated pro- 
duction 



fl fl/zmr /"tmax r\(T i 

a = dr dy dt £ f a/A (x aj M 2 a ) f b/B (x b , M 2 b ) (15) 

Jm 2 /S J-l/21nr Jt min ™ 

can be obtained by convolving the partonic cross sections da a b/dt presented in Sec. Ill Al 
with the parton density functions (PDFs) f a ,b/A,B at the factorization scale M a ^. Since 
the PDFs vanish rapidly, as the longitudinal momentum fractions x a j> of the partons a, b 
in the external hadrons A, B approach unity, the available partonic center-of-mass energy 
s = x a XbS represents only a fraction r = x a Xb of hadronic center-of-mass energy S, and the 
experimentally accessible mass range for searches of new SUSY particles is naturally limited. 
We consider the initial gluons and five light quarks to be massless and denote the average 
final state mass by m. At the LHC, both A and B represent protons, which will collide with 
\/~S = 14 TeV starting in 2008, whereas at the Tevatron, B represents an anti-proton beam 
with vS = 1.8 TeV at the completed Run I and 1.96 TeV at the current Run II. 

Numerical predictions for single-gravitino hadroproduction cross sections at Run I of 

n 

the Tevatron have been presented in Figs. 4, 5, and 6 of [12] as a function of m~ g and 
mq, respectively. We have verified these results by fixing the strong coupling to its world 
average value a s (Mz) = 0.118 and convolving our partonic cross sections in Sec. Ill Al with 
the (nowadays obsolete) set of PDFs of jzj, evolved from the starting scale Q = 2 GeV 
and using a value of A^o 5 = 144 MeV to the factorization scale M n = Mb = rrig^. In 
the following, we will, however, use the modern PDFs of CTEQ6L1 [23], which correspond 
to a one-loop running of the strong coupling a s (/i) = ^/(47r) and a QCD scale parameter 



of Alq 5 = 165 MeV, derived from the world-average value of a s (Mz) = 0.118 19j. The 



renormalization scale /i and the factorization scales M a ^ will be fixed to the average particle 
mass m = (m^ + mq^)/2 in the final state. 

In Fig. [HI the total cross section of the associated production of gravitinos and gluinos 
or squarks is shown for the GMSB benchmark scenarios SPS 7 (top) and 8 (bottom) and a 
gravitino mass of rriQ = 10 -5 eV as a function of the effective SUSY-breaking scale A. While 
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FIG. 8: Total cross sections of gravitino and gluino/squark associated production at Run II of the 
Tevatron for = 10 -5 eV and the GMSB benchmark slopes SPS 7 (top) and SPS 8 (bottom) as 
a function of the effective SUSY-breaking scale A. 
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the quark-gluon initiated production of gravitinos and squarks can contribute significantly 
to the total event sample for SPS 7, where m g - < nig, the largest contribution at SPS 7, 
and even more so at SPS 8, where nig < Trig, comes from the subprocess qq — > Gg. This is 
of course due to the large quark- ant iquark luminosity at the Tevatron. At Run II, where 
the integrated luminosity has already reached 2.6 fb -1 and is expected to increase to 4.4 

the CDF and DO experiments should be 



to 8.8 fb" 1 until the final shutdown in 2009 |l_ 
able to discover light gravitinos with masses up to 10 -5 eV in events with a single jet and 
large missing transverse energy for values of A above the current exclusion limits (25 and 70 
TeV for SPS 7 and 8, respectively). If we assume the large- Et monojet signal to have very 
little SM background (see Fig. [9]below) and thus to be experimentally identifiable with high 
efficiency, we can define the visible region by the point where the total cross section falls to 
1 fb, so that only a few events will be recorded. The discovery reach then extends up to the 
benchmark point (40 TeV) at SPS 7, where nig = 950 GeV and Trig = 890 GeV, and even 
up to 120 TeV at SPS 8, where m g = 1000 GeV and Trig = 1300 GeV. For lighter gravitino 
masses, the total cross section scales trivially according to Eqs. ([3]), (J5j), and (J7J), i.e. with 
the inverse of the squared gravitino mass. 

For a detailed, model-independent account of the experimentally identifiable monojet 
signal of the associated production of a gravitino with a squark or gluino, we have to include 
the decay of the latter into a second gravitino and an observed jet. If we continue to neglect 
rriQ, except in the coupling, the cross section for a massless three-body massless final state 
can be written as 

da = ^//aM(^,M a 2 )dx a / b/B (x 6 ,M 6 2 )dx fe pWT^ 3 (27r) 4 5 4 L+J*-J>) A (2 ffi^ 

= ^//aM(^,M a 2 )/ fc/B (a; b ,M 2 )^ (16) 

where px 1 represents the observed jet transverse momentum, which is balanced by the missing 
transverse momentum of the two gravitinos. Since the squark or gluino decay width is of the 
order of he / 1 fm = 0.2 GeV (see Sec. IIIIBj) . we can apply the narrow- width approximation 
to rewrite the squared and averaged 2 — 3 scattering matrix element as 

2 

I »/f|2 

+2 



\m\Ls = imiL 2 



i 



S i2 - m\g + irXlg-gYg^g 



\M\l 



mL 2 nS{Sl2 - m ^ \M\u. d7) 

m q,g 1 q,g 
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By fixing the azimuthal angle of the observed jet to 0i = 0, the squared invariant mass of 
the intermediate squark/gluino propagator becomes 



S12 = 2pt 1 Vt 2 (cosh 771 cosh 772 — sinh 771 sinh 77 2 — cos fa), (18) 

so that 



8 



5(s 12 - m~ 



2 — arccos cosh 771 cosh 772 — sinh 771 sinh 772 



?.) = 

p Tl PT 2 \ sin 02 1 



m? - 

9-9 
ZPT-lPT.2 



(19) 



The three-body cross section 

da = [ x a f a/A (x a , M^)x b f b/B (x b , M%)^-BR(q, g -> GAT) 2 dp Tl d77idpT 2 d77 2 d773 (20) 

J Qt 7T Sill ©2 



can then be expressed in terms of the squared and averaged 2 — > 2 production cross section of 
gravitinos and squarks or gluinos, da/dt (see Sec. Ill Aj) . and the squark or gluino branching 



ratio BR(g, g — > GX) into a gravitino and a jet (see Sec. IIIBI) . 

The transverse-momentum spectrum of the observed jet, which is equivalent to the miss- 
ing transverse-momentum spectrum, is shown in Fig. [9] for the same gravitino mass of 10~ 5 
eV as in Fig. [HJ At this point, squarks and gluinos always decay into gravitinos and jets (see 
Fig. [6]). Their decay widths vary between 4 and 17% of their masses, as these increase from 
the current exclusion limit of 250 GeV to 350 GeV, so that the narrow-width approximation 
is always justified. All three spectra peak at values slightly below half of the squark/gluino 
mass, as expected from kinematic considerations. The main SM background, which comes 
from the associated production of a jet and a Z-boson, followed by an invisible Z-decay, 
peaks at roughly half the Z-boson mass. It can be eliminated by cutting on the invisible (or 
jet) transverse momentum at values around 100 GeV. 

In a recent analysis of Run I Tevatron data, the CDF collaboration have examined events 
with a single jet and a missing transverse energy of at least 100 GeV 25j- From Fig. [9] it is 
clear that, while this cut eliminates basically all of the soft-QCD and other Standard Model 
backgrounds, only little signal cross section is lost. For an optimized missing transverse 
energy cut of 175 GeV, the CDF collaboration found an upper limit of the gravitino cross 
section of 3.1 pb, corresponding to a gravitino mass of at most 1.1 ■ 10 -5 eV. Note however 
that this analysis was done under the assumption that all other supersymmetric particles 



are heavy 26]. 
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FIG. 9: Transverse-momentum spectra of the observed jet in gravitino production at Run II of the 
Tevatron for tUq = 10 -5 eV and three different squark/gluino masses. Also shown is the main SM 
background from invisible Z-boson decays. 



We therefore repeat the CDF analysis for general SUSY scenarios, using the 95% con- 
fidence limits on the product of the acceptance (A) times the signal cross section (a) as a 
function of missing- .&r as published in Fig. 3 of Ref. [25j. These limits are divided by the 
detector acceptance (A = 0.4) for the selected data sample, assumed to have little depen- 
dence on the missing- j^]. The dependence of the number of simulated signal events on 
the varying missing- Ex cut is then taken into account explicitly by integrating Eq. (I2"0"j) over 
Vtx — 100 ... 300 GeV. We also implement the experimental requirement that at least one 
jet lies in the central region, \rj\ < 0.7, since this cut is stricter than the additional CDF-cut 
on the hardest jet to lie in \r)\ < 2.4, and our parton-level analysis has only one jet. We have 
verified that the number of events with a hard jet in the region 0.7 < \tj\ < 2.4 is indeed 
negligible. 



From our confirmation of Figs. 5 and 6 in 



12] we know that the two production processes 



that involve both the squark and gluino mass (qq — > Gg and qg Gq) are bounded from 
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FIG. 10: Exclusion contour (full curve) in the m^, — nig^ plane derived from the CDF acceptance 
times cross section limits for events with a single jet and varying missing transverse energy [3, H]- 
Also shown are the validity region of the narrow-width approximation (NWA, above the dotted 
curve), the CDF limit obtained for very heavy squarks and gluinos (dashed line), and the unitarity 
limit (dot-dashed line). 



below for nig = nig, while the third production process gg — > Gg depends only on the gluino 
mass. We sum therefore over all three production subprocesses with rrig = m g and take BR 
= 1 for ttlq < 10~ 4 eV according to Fig. [61 

By imposing always the strongest experimental experimental limit on A a ($t) of Fig. 3 



in Ref. 



251 ]. divided by A = 0.4 [27|, on the correspondingly integrated cross section, Eq. 



( I20|) . we obtain a contour in the — m g ^ plane, which is shown in Fig. [10] (full curve). 
For light squark and gluino masses of 200 GeV, we find a gravitino mass limit of 2 • 10 -5 
eV that is very similar to that found by CDF for very heavy squarks and gluinos (dashed 
line). The limit in Fig. [T0l degrades slowly to 4- 10 -6 eV as the squark/gluino mass increases 
to 700 GeV, i.e. as it approaches the center-of-mass energy available at the Tevatron and 
the theoretical cross section falls. At the same time, the squark/gluino width increases 
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and eventually passes the value of 1/4 of the squark/gluino mass (dotted curve in Fig. [TUI) . 
Our results, obtained in the narrow-width approximation for single-gravitino production 
in association with relatively light squarks and gluinos, are thus complementary to those 
obtained by CDF for very heavy squarks and gluinos and double-gravitino production |25|. 
We have also checked that for our analysis, which is based on a single-gravitino effective 
Lagrangian, tree-level unitarity is always satisfied [28|, since 

G > 0.3 9 - — (21) 

- inn n^r V / 



10- 6 eV - ' 100 GeV 

for a critical energy corresponding to the Tevatron center-of-mass energy (dot-dashed curve 
in Fig. HO]). 

Two other experimental analyses of monojet signals at the Tevatron have been published, 
one based on 78.8 pb _1 of Run-I data by the DO collaboration [29] and one based on 368 
pb _1 of Run-II data by the CDF-collaboration 3(3]. However, both analyses are interpreted 
with extra-dimensional models and directly present limits on the number of these extra 
dimensions and on the corresponding fundamental Planck scale. It would be interesting to 
re-interpret these analyses in the context of gravitino production. The CDF analysis quotes 
indeed a model- independent limit on signal events (or signal cross section times acceptance), 
but neither publication gives numerical values for detector acceptances. These were also not 
available from the collaborations upon request, so that we can at this point not deduce 
independent gravitino mass limits from these data. 



D. LHC 

The high center-of-mass energy of v^S* = 14 TeV and the large luminosity of initially 
10 fb^ 1 and finally 300 fb _1 available at the LHC will provide the opportunity to test the 
soft SUSY-breaking hypothesis up to the multi-TeV range. This general remark remains 
true for the gauge-mediated SUSY-breaking scenarios SPS 7 (top) and 8 (bottom) with a 
gravitino LSP that we consider in Fig. [TTJ In this figure, we show the total cross sections 
of gravitino and gluino/squark associated production at the LHC for itlq = 10~ 4 eV and 
the three different partonic gravitino production processes discussed above. Their hierarchy 
is now opposite to the one at the Tevatron, i.e. it is the gluon luminosity that dominates 
and no longer the quark-antiquark luminosity. At SPS 7, where m ? - < m§, squarks are 
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FIG. 11: Total cross sections of gravitino and gluino/squark associated production at the LHC for 
ttiq = 10 -4 eV and the GMSB benchmark slopes SPS 7 (top) and SPS 8 (bottom) as a function 
of the effective SUSY-breaking scale A. 
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FIG. 12: Transverse-momentum spectra of the observed jet in gravitino production at the LHC 
for SPS 7 (full line) and SPS 8 (dashed line). 

produced more copiously than gluinos, wheres the inverse is true at SPS 8. In both cases, 
the discovery reach extends to values of the effective SUSY-breaking scale A far above 
the actual benchmark points. It is clear that a striking monojet signal with large missing 
transverse energy could be discovered rapidly after the start-up of the LHC and with rather 
low luminosity. The only existing previous analysis for gravitino production at the LHC, 
based on the dominating subprocess gg — > Gg only, assumed a slightly higher LHC center- 
of-mass energy of y/S =16 TeV and fixed squark and gluino masses of 2 TeV and 750 GeV, 
respectively [9|. The authors computed a monojet cross section of similar size as the one 
shown in the lower part of Fig. [TT] (8 pb for = 10~ 4 eV) for comparable squark and 
gluino masses. 

The transverse-momentum spectra at the LHC are shown in Fig. [T2] for the two GMSB 
benchmark points and a gravitino mass of uiq = 10 -4 eV. While the spectra peak again 
at values slightly below half of the squark/gluino masses, as was already the case at the 
Tevatron (see Fig. [9]), they extend now to much larger values of px — 1200 GeV. Note also 
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that the absolute Tevatron cross section was only of similar size since it had been calculated 
with a smaller gravitino mass of m<§ = 10 -6 eV. In the high-luminosity phase at the LHC, 
the missing transverse-energy signal will be degraded by pile-up events. In the ATLAS 
experiment, pile-up can be eliminated and acceptable trigger rates in the kHz-range can be 
obtained for missing- E? thresholds of 60/120 GeV at low/high luminosity. These values can 
be two times lower if an addition hard jet of Et > 100 GeV is required 31[. As can be seen 



from Fig. [121 the signal cross section will be affected very little by these cuts. 



IV. CONCLUSION 



In this paper, we have derived the Feynman rules for light gravitino production and 
decay from an effective supergravity Lagrangian in four-component, non-derivative form and 
computed analytical partonic cross sections and branching ratios involving interactions of 
gravitinos, gluinos, and squarks. Special emphasis has been put on the gauge-independence 
of the results, the contributions of quartic vertices, and a comparison with results obtained 
previously with a derivative Lagrangian and those obtained in the high-energy limit. 

Using the narrow-width approximation, we combined the associated gravitino- 
squark/gluino production cross sections with the subsequent decay of the squarks/gluinos 
into quarks/gluons and a second gravitino. This enabled us to perform extensive numerical 
studies of branching ratios, of the total hadronic cross sections at the Tevatron and the 
LHC, and of the corresponding transverse-energy spectra of the single observed jet, then 
to impose experimental cuts on the latter, and finally to derive a new and robust exclusion 
contour in the m^/m^-plane from the latest CDF mono jet cross section limit. 

Our Tevatron exclusion contour implies that gravitinos with masses below 2 ■ 10 -5 to 
1 • 10~ 5 eV are excluded for squark/gluino-masses below 200 and 500 GeV, respectively. 
These limits are complementary to the one obtained by the CDF collaboration, 1.1 • 10~ 5 
eV, obtained under the assumption of very heavy squarks and gluinos. 

For the LHC, we conclude that SUSY scenarios with light gravitinos, such as the GMSB 
benchmark slopes SPS 7 and 8, will lead to a striking monojet signal very quickly after its 
startup in 2008 and already with low luminosity. The missing- Et and jet trigger thresholds 
foreseen by the ATLAS collaboration are perfectly suitable also in these scenarios for an 
efficient background reduction without affecting the signal in a significant way. 
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APPENDIX A: EFFECTIVE LAGRANGIAN FOR SINGLE GOLDSTINOS 



We start from the effec 



tion of a light gravitino 



;ive Lagrangian in four- component notation for the single interac 



IB. 



32. 



whose longitudinal (spin-1/2) goldstino components 
can couple to the matter and gauge supermultiplets with enhanced (electroweak) strength 
lL The corres ponding effective Lagrangian in two-component notation can be found in 



14 



34 



35 



36 



37 



38 



39 



401 ]. while interactions involving several external goldstinos or 



goldstino propagators and the scalar /pseudo-scalar superpartners of the goldstinos, the so- 
called sgoldstinos, have been derived in four-component notation, e.g., in 33l. 141]. 

By correctly translating the effective Lagrangian for light gravitinos in two-component 



notation 



n 



391 ] into the four- component notation of the traditional SUSY-QCD Lagrangian 
we obtain the following non-derivative couplings for the interactions of Majorana- 
fermionic goldstinos ip with Dirac-fermionic quarks \i complex-scalar squarks 0, vector- 
bosonic gluons A^, and Majorana-fermionic gluinos A a : 



+ 



irrtn 



g s m~ g 



V6Mm, 



(Al) 



Here, M = (%-kG n )- 1 / 2 = 2.435 • 10 18 GeV is the reduced Planck mass, rriQ is the gravitino 
mass, which is related to the supersymmetry breaking vacuum expectation value (F) in 
canonical normalization by vtlq = (F)/(v / 3M), and g s is the strong gauge coupling. TP- 
axe the generators of the SU(Nc = 3) color symmetry group with antisymmetric structure 
constants f abc and Casimir operator Cf = 4/3, and Pl.r — (1 -F Ts)/2 are the chirality 
projection operators. The squark and gluino masses will be denoted trig and nig. Since the 
top quark density in hadrons is small, we can neglect the masses m q of the five light quarks 
at high collision energies and consider the corresponding left- and right-handed squarks to 
be mass-degenerate. 
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The effective theory contains the same couplings as the full theory 43|, l44j . except for 



the quark-gluon-gravitino-squark vertex, which would violate the gauge symmetry in the 



effective theory 14| , but was erroneously kept in the alternate Feynman rules of 



m. These 



were, however, not used there to calculate cross sections. In contrast, there is a new four- 



particle vertex [jj], |39| , the gravitino-gluino-squark-squark vertex, which has been overlooked 
in all other cited references, but is neither relevant for our analysis nor for the one in 12]. 
Attention must also be paid to the assignment of factors of i in Eq. flAll) . if the interference 



terms in 



121 ] between scalar and gauge boson exchanges are to be correctly reproduced. 



In the effective Lagrangian, all vertices are proportional to SUSY-breaking mass terms, 



i.e. m~ — ml and m g . In particular, the Yukawa coupling of the goldstino can be obtained 



m- 



from that of the gluino by the replacement 

a T a - - 

while the goldstino-gluon-gluino coupling can be obtained by the replacement 



g s f ahc r - i 



m 9 



2V6M- 



ab 



(A2) 



(A3) 



G 



with P representing the incoming gluon-momentum. At high energies, contributions involv- 
ing the cubic goldstino-quark-squark coupling are suppressed relative to the gluino contri- 
bution due to the higher mass-dimension of the coupling. 



APPENDIX B: FEYNMAN RULES FOR SINGLE GOLDSTINOS 

In order to derive the Feynman rules needed for the hadronic production and decay of 
single goldstinos, we multiply the effective Lagrangian in Eq. ( 1A1I) with a factor of i, perform 
a Fourier transformation, and take the functional derivative with respect to the external 
fields. Denoting the incoming four-momentum by P, the chirality projection operators 
by Pl,r = (It 7s)/2, Lorentz indices by /x, u, and color indices of the fundamental 
(adjoint) representation of the color symmetry group SU(3) by i, j, ... (a, b, ...), we obtain 
the following interaction vertices: 




P,a,n 



+ i 



m-g 



sab 



2V6Mm r 



.7" 



(Bl) 
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a, ii 



+ 



2 v / 6Mm (5 



9s f abc [7 M , Tl 



(B2) 



2 2 

m~ — ra„ 
=F« — p 2.<J«P; 



ij *L,R 



(B3) 



■- J 



9 % Pr,l 



(B4) 




' 9s T { 



(B5) 



Here, the arrows on (s)quark lines indicate flavor flow, while the Majorana nature of 
gravitinos and gluinos requires the fermion flow to be fixed arbitrarily 45j. These Feynman 



rules have been implemented into the computer algebra program FeynArts 46|, |47|], and the 



corresponding model file is available from the authors upon request. The derivative forms of 
the (s)goldstino interaction vertices in two-component form have already been implemented 
some time ago into the program CompHEP [48 ] . 



Our Feynman rules differ from those presented in appendix A. 3. 3 of 49(, which have 
also been derived from [141 ] . by a factor of % in the first two vertices, apparently due to a 
misinterpretation of the two-component tensor 50]. The Feynman rules in 5J| differ 
from ours in addition by an irrelevant global factor of i. The usual SUSY-QCD vertices and 
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propagators can be found, e.g., in 
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